A theoretical model for attachment lifetimes of kinetochore-microtubules: Mechano-kinetic 

"catch-bond" mechanism for error-correction. 
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Before cell division, two identical copies of chromosomes are pulled apart by microtubule (MT) filaments 
that approach the chromosomes from the opposite poles a mitotic spindle. Connection between the MTs and 
the chromosomes are mediated by a molecular complex called kinetochore. An externally applied tension can 
lead to detachment of the MTs from the kinetochore; the mean lifetime of such an attachment is essentially a 
mean first-passage time. In their in-vitro pioneering single-kinetochore experiments, Akiyoshi et al. (Nature 
468, 576 (2010)), observed that the mean lifetimes of reconstituted MT-kinetochore attachments vary non- 
monotonically with increasing tension. The counter-intuitive stabilization of the attachments by small load 
forces was interpreted in terms of a catch-bond- like mechanism based on a phenomenological 2- state kinetic 
model. Here we develop the first detailed microscopic model for studying the dependence of the lifetime of 
the MT-kinetochore attachment on (a) the structure, (b) energetics, and (c) kinetics of the coupling. The catch- 
bond-like mechanism emerges naturally from this model. Moreover, in-silico experiments on this model reveal 
further interesting phenomena, arising from the subtle effects of competing sub-processes, which are likely to 
motivate new experiments in this emerging area of single-particle biophysics. 



Chromosome segregation by the mitotic spindle is one 
of the most important intracellular processes in eukaryotic 
cells Ill2r4l5l1 . Connections between chromosomes and micro- 
tubules (MT) are mediated by kinetochores, which are com- 
plex macromolecular structures H16l - l20ll . Due to the difficul- 
ties of isolating kinetochores from cells, the identification and 
spatial organization of the molecular components of kineto- 
chores has posed significant challenges. Recent high reso- 
lution imaging has provided an indication of the distribution 
of these components and even their stoichiometries 1321 - 13 411 . 
Kinetochores form dynamic, and yet sufficiently strong, cou- 
pling with MTs that undergo stochastic transitions between 
growth and shortening. This interaction between kinetochore 
elements and the attached kinetochore MT (kMT) generates 
movement of the chromosome. While a possible architec- 
ture of this nano-device is beginning to emerge, the mech- 
anism by which it couples forces from kMT polymeriza- 
tion/depolymerization with chromosome movement is a major 
unresolved question with significant implications 
A fundamental biophysical question in this context is: how the 
dynamics of the kMTs and externally applied tension (load 
force) affect the stability of the kMT-kinetochore coupling. 

Recent in-vitro experiments with reconstituted kinetochore- 
MT attachments in budding yeast @, [^TJ] have provided evi- 
dence that MT kinetics and load forces can combine in un- 
expected ways. Strikingly, it was found that a limited range 
of forces can be more favorable for maintaing kinetochore at- 
tachment, whereby load selectively stabilizes attachment 112111 . 

In this letter, we develop a detailed microscopic model that, 
to our knowledge, is the first theoretical analysis of this phe- 
nomenon, at the single kinetochore level. The mean lifetime 
of the MT-kinetochore attachment is essentially a mean first- 
passage time Hh . Calculating this mean first-passage time us- 
ing our microscopic model, we investigate the dependence of 
the mean attachment lifetime on (i) the structure, (ii) energet- 
ics, and (iii) kinetics of the MT-kinetochore coupler. Akiyoshi 



et al. L21] argued that their counter-intuitive data are "rem- 
iniscent of 'catch-bonds'" that can be explained in terms of 
a phenomenological two- state kinetic model. A catch-bond- 
like mechanism emerges naturally in the theoretical frame- 
work of our microscopic model as a consequence of the force- 
sensitivity of kMT depolymerization. Our results reveal wider 
varieties of trends of variation of the attachment time than 
those observed by Akiyoshi et al. We also indicate pos- 
sible adaptations of the experimental techniques of Akiyoshi 
et al. lElh that may be appropriate for testing our new predic- 
tions. 

Almost all the theoretical models of MT-kinetochore cou- 
pling i23r - [27il29ll are based exclusively on one of the two ma- 
jor mechanisms for force generation. In the biased-diffusion 
model, initially proposed by Hill [23], the plus end of a kMT 
is assumed to be surrounded by a coaxial "sleeve" the inner 
surface of which is composed of several binding elements 
that bind specific kMT sites. The one-dimensional Brown- 
ian motion of the sleeve along the axis of the kMT is bi- 
ased to increase overlap, because a larger number of kMT- 
sleeve bindings lowers the total energy of the system. The 
interplay of this biased diffusion and the depolymerization 
of the kMT gives rise to the pull exerted by the coupler on 
the kinetochore. An alternative coupling mechanism is based 
on the "power stroke" exerted on a rigid rin g by the curling 



protofilament tips of a depolymerizing MT [261 |27|] . There 
is increasing recent structural evidence that kinetochores in- 
deed engage kMTs through multivalent attachments that move 
along microtubules • Therefore, a biased diffusion re- 

mains a valid candidate mechanism for MT-kinetochore cou- 
pling. But this evidences does not necessarily exclude a role 
of the well known curled tips of depolymerising MTs in the 
MT-kinetochore coupling. 

In contrast to most of the earlier theoretical work, the model 
we propose here is "unified" in the sense that it incorporates 
the key features of both these types of models. In our model 
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the main elements of the biased diffusion model are treated 
explicitly. Moreover, the curling of the MT protofilament tips, 
a key feature of the power-stroke model, is captured implic- 
itly by assuming a tension-induced slowing down of depoly- 
merisation which is known to arise from the tension-induced 
suppression of the curling. Although force-mediated kMT al- 
teration has been discussed in the literature modeling 
of these effects in the context of biased diffusion has been lim- 
ited. 

The theory we develop here is an extension of a one- 
dimensional force-based model of a kMT-kinetochore "cou- 
pler" [29]. In this extended version, the coupler is composed 
of multiple passive kinetochore elements that bind kMTs via 
a generalized biased-diffusion mechanism, along with active 
kinetochore force generators, all of which maintain dynamic 
attachment with shortening/growing attached microtubules. 

We begin with the simplified special version of our model 
that includes only the passive binders; later in this letter we 
incorporate also the active force generators. In the first simpli- 
fied version, the kinetochore coupler is modelled as a collec- 
tion of binder element heads which represent the core binding 
area of a kinetochore. 

The length of the overlap between the kMT and the cou- 
pler is denoted by x (see Fig [I}. Increasing the overlap be- 
tween binders and the lattice is energetically favorable. As 
in previous work fl29[] , we assume that each binder head en- 
gages with the kMT by obeying a unit energy function (j>b{x), 
which has two key parameters: a measures free energy drop 
due to binder affinity for the kMT lattice, and b describes the 
activation barrier for transitions between specific kMT lattice 
binding sites. In other words, a is a measure of the strength 
of the kMT-binder affinity while b is a measure of the "rough- 
ness" of the kMT-coupler interface. The total potential energy 
function is given by 



Binding Energy/Binder 



^b(x) = ^2(j) b (x - ns) 



(1) 



where s is the spacing between consecutive coupler binders 
(see Fig [I]). Binder spacing is an arbitrary parameter (SI). 
Here we set s = I, where £ is the distance between consecu- 
tive kMT binding sites. 

The coupler overlap velocity, dx/dt is then given by the 
following stochastic differential equation 



dx(t) 



WFdt 



j [-%(x) - F load ] dt 



(2) 



• ldN r {t) + y/2k B T/tdW(t), 



where the constant Fi oa( j is the external opposing load force 
on the coupler and £ is the effective drag coefficient. dW(t) 
accounts for the thermal diffusion of the coupler on the lattice, 
and N(t) is poisson counting processes describing the kMT 
dynamics with intensity rate r. 

An opposing tension tends to decrease the overlap between 
the kMT and the coupler. We assume that the coupler under 
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FIG. 1. Diagram of model components, x = is the kMT entry- 
point and x = L is the maximum overlap. Mechano-kinetic mod- 
ification leads to amplification of energy barriers for binding in a 
'catch-bond' type mechanism. 



tension also suppresses the curvature of the tips of the kMT 
protofilaments so that the splaying tips of the kMT become 
confined within the coupler (we visualize the coupler protein 
mesh work as a children's finger trap toy, see Fig [T}. Such 
force-induced suppression of the curvature of the depolymer- 
izing protofilament tips, in turn, causes reduction of the rate 
of kMT depolymerization. This proposed scenario is consis- 
tent with the experimental observations of [|2lll . Thus, kMT 
depolymerization rate is assumed to be a decreasing function 
of the load tension and the functional form of this dependence 
is assumed to be 



-A Fi oa d 



(3) 



where the parameter A characterizes the extent of the effect of 
a given load tension on f3. The kMT polymerization rate a, 
however, is assumed to be independent of load. The rate func- 
tions are chosen in order to allow for the rate r = a — /3(F) 
of the kMT tip dynamics to transition from a catastrophe to 
a rescue state in a load-dependent manner in agreement with 
observations in [21] (SI). We define the breaking load Fb rea k to 
be the strength of the tension for which the mean attachment 
time is less than 1 min; the qualitative conclusions drawn on 
the basis of this definition do not depend sensitively on this 
choice. 

The lifetime of a MT-kinetochore attachment is defined 
here mathematically as the time taken by the kMT, that is ini- 
tially at x = L, to reach x = (the coupler entry point) for the 
first time. Since this time fluctuates from one MT-kinetochore 
attachment to another, we calculate the mean lifetime. The 
mean attachment lifetime is essentially a mean first passage 
time [1], which we calculate using standard methods (SI). 
By a combination of analytical and numerical techniques, we 
study the trends of variation of the mean lifetime with (a) the 
strength of the externally applied load force, as well as, the 
(b) microscopic structure, (c) energetics, and (d) kinetics of 
the coupler. iV& is characteristic of the structure of the coupler 
(coupler length) whereas its energetics depend on ^5 (i.e., on 
the parameters a, b) and Fi oad ; the stochastic kinetics are in- 
fluenced by the interplay of forces arising from the potential 
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landscape, random Brownian forces, and by the kMT poly- 
merization / depolymerization kinetics. 

It is difficult to derive an exact analytical expression for 
the mean first passage time related to eq. (0). Therefore, we 
explore two limiting cases for which explicit approximate so- 
lutions can be obtained: (a) Slippery regime (i.e., low-friction 
regime) where b « ksT; in this regime the coupler can eas- 
ily rearrange its position relative to the kMT, (b) Strong fric- 
tion regime where b » /c#T; in this regime diffusion inside 
the binder is practically non-existent and MT growth/ shorten- 
ing rates are large compared to all other processes. Stronger 
friction weakens the ability of the coupler to quickly adjust its 
position with the variation of the length of a dynamic MT. For 
sufficiently large b, the coupler becomes static ll29ll . 

In the slippery regime, the mean lifetime is (see SI for the 
derivation) 



T(L) 



L 2 exp (—w) — 1 



D 



(4) 



where w = L(-a/£ + F load + £f3 max exp(-\F load )£)/k B T is 
a dimensionless work quantity. 

In the strong friction regime, the mean lifetime is (SI) 



T(L) 



L 



-A Fi 0; 



(5) 



In Fig [2 we show plots of mean lifetimes obtained by com- 
puter simulations for various parameter values, in addition to 
the expressions © and ©. Over a significant regime of phys- 
ically relevant parameter values, our model gives rise to a n on 
monotonic variation of the mean lifetime with the load ten- 
sion; this is consistent with the experimental observation of 
Akiyoshi et al. lEHl . We have also explored other parame- 
ter regimes to understand kMT-kinetochore detachment phe- 
nomenon in further detail. The results indicate the possibility 
of other distinct trends of variation of the mean lifetime with 
tension that might be detectable in experiments under condi- 
tions different from those used by Akiyoshi et al. 
Results for the slippery regime: FigsEl and[3]\ show that in 
the slippery regime the trend of variation of the mean life 
time on the load tension depends sensitively on the binding 
energy (provided by a or A5). At sufficiently low binding 
energies, attachment times can be sensitive to depolymeriza- 
tion. Under these conditions, if depolymerization slows down 
by the load before the coupler breaks, then attachment times 
are non-monotonic and essentially follow /3(F) (A5 = 32 in 
Fig. [2]). The peak value of the mean attachment time depends 
on A = 1/F C . Thus, for observing the non monotonic varia- 
tion of the mean life time, A must be chosen so as to satisfy 
the requirement F c < Fb rea k. Consequently, for F- indepen- 
dent f3 (which corresponds to the special case A = 0), the 
mean lifetimes decrease monotonically with increasing load 
tension, in agreement with 112 ill (SI). As our results establish, 
the range of breaking load Fb rea k can be easily adjusted for 
these model couplers by increasing the binding affinity for the 
lattice using a or (two cases shown in Fig. 0. However, 
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FIG. 2. Mean attachment time versus load force. Blue, orange 
and red correspond to the slippery, intermediate and strong friction 
regimes, respectively. Solid blue line is obtained from eq. J4j with 
A = 3pN~\a = 40,Ana X = 120, N b = 65, a = 0Ak B T,b = 
0.001a; dashed blue line is obtained numerically for the same pa- 
rameters but with Nb — 32. Parameters for solid orange triangles 
a = 50, /3 m ax = 350, N b = 45, a = 0.5/c s T, b = 0.04a; those for 
solid orange squares are same except that a = 0.6/csT, b = 0.2a. 
Solid red curve is obtained from eq. (|5} with (3 mSLX = 350, a = 
0,A = 0.4. 



this enhanced stability comes at a price: in the stable slippery 
regime (A5 = 65 in Fig. 0, the lifetimes for low loads in- 
crease beyond observable ranges and sharply decrease close 
to the breaking loads. 

Results for the strong friction regime: In this regime the dif- 
fusive motion of the coupler is made practically impossible 
by the condition b <C ksT. Moreover, increasing load ten- 
sion cause stronger suppression of the kMT depolymerisation. 
Consequently, in this regime, the mean life time increases 
monotonically with increasing load tension (see Figj2j). 
Results for the intermediate regime: 

We investigated the intermediate regimes numerically. The 
data in Fig. [3] establish that, in this regime, molecular friction 
makes it harder for the kMT to exit the coupler under load, 
while it also further enhances the dependence of attachment 
time on kMT depolymerization. A delicate balance between 
these effects strong friction and force leads to the observed 
trend of variation of mean life time with load force in this 
regime. 

The attachment regimes that result from our model reveal 
competing effects at the kinetochore sites. On one hand, load 
forces provide a pull that can detach the kMT from the kineto- 
chore coupler. On the other hand, load force slows down de- 
polymerization of the kMT, thereby slowing down the exit of 
its tip from the coupler which, effectively, counters the pulling 
effects of the load. In the conflict between these two opposing 
effects, as indicated by our data in the Figs. [2]and[3j the inter- 
nal coupler friction might be the ultimate determinant of the 
emergent behavior. In low and intermediate friction ranges, 
we distinguish a clear range of forces for which coupling is 
selectively favored. This finding supports a force-mediated 



4 



^a=0.45k B T,b=0.04a 
^a=0.5k B T,b=0.04a 
^a=0.6k T,b=0.2a 




A ^.Passive Binder A Minus-End Motor 



Force (pN) 



Force (pN) 



FIG. 3. Mean attachment times for various coupler parameter 
regimes. A. Mean attachment times in the slippery regime with dif- 
ferent kMT kinetic rates for a short coupler with Nb — 32. Common 
parameters are A = 3 pN -1 , a = 0.4 k s T, k = 0.001a. B. Mean 
attachment time variations for an intermediate regime coupler with 
Nb = 45 binders. An increase in the binding energies a causes the 
tails of the mean attachment times to expand in larger force regions. 
Common parameters are A = 3 pN -1 , (3max — 350 s _1 , a — 50 
s _1 . Error bars mark standard deviation. 



selective stabilization of kMT at the kinetochore coupler. We 
find that non monotonic variation of the mean life time with 
the tension, which resembles a a catch-bond mechanism sim- 
ilar to[21] is just one of the possible responses of the kMT- 
kinetochore coupler. It might be possible to observe the other 
types of theoretically predicted responses by creating the cor- 
responding required conditions in the in-vitro experiments. 
Such conditions may be facilitated, for example, by biochem- 
ical modifications of the Ndc80 complexes at kinetochores 
which are close candidates for our multivalent pas- 
sive binders. 

Effects of force -generating motor proteins in the coupler 
In addition to the passive kMT binders, active force gener- 
ating components also play an important role in maintaining 
and regulating kMT-kinetochore coupling. Among the active 
force generators, cytoskeletal motorproteins are believed to 
make the dominant contribution IeI-S [ZD • Using recent struc- 
tural data Jsl Ei] we create a hybrid coupler, where the out- 
ermost layer is composed of passive components and the in- 
nermost layer closest to the chromosome is composed of an 
active interface (see Fig[6}\). 

To include these active components, we add an active force 
term in our model 



F A (x) = d m (x)(n-f- - n+f+) 



(6) 



where d m (x) measures the x-dependent length of the active 
kinetochore interface. The parameters n_, n + denote the av- 
erage number of minus and plus end motors per unit length of 
MT embedded in the kinetochore structure while /_ and /+ 
denote the force generated by a single minus-end and plus-end 
directed motor, respectively. For each motor we postulate a 
linear force- velocity relation in agreement with previous work 

A key feature of the hybrid coupler is that besides modify- 
ing the force balance of the coupler, active components also 
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FIG. 4. A. Diagram of the architecture of a hybrid coupler. B. Nu- 
merical simulation results of average speed versus force for a stably 
bound coupler with varying densities of minus-end directed motors. 
Inset. Coupler stall forces for each motor density. C. Average speed 
versus force for a stably bound coupler with varying densities of plus- 
end directed motors. D. Breaking load calculations for the coupler 
shown in panel B-C for two binding energy values. The parameters 
used are Nb = 52, A = 3, f3 ma x = 130 s _1 , a = 20 s _1 , a = kgT, 
k 0.01a. 



affect the internal friction coefficient (see SI for details). 

In Fig. [6] we plot force- velocity relations and F break for a 
hybrid coupler with varying numbers of motors at the active 
interface. The non-monotonic nature of mean attachment is 
preserved when the active components are added. When the 
motors oppose load force, stability of the kMT-kinetochore 
coupling is enhanced, with dynein (or protofilament curling) 
being a strong candidate for this role. This effect is supported 
by the increased coupler breaking loads as n_ increases, Fig. 
[6t). The active interface stabilizing effect of minus-end mo- 
tors is particularly amplified when the binding energy of pas- 
sive components is weakened, Fig [(J). Despite the stability of 
the attachment being improved, higher numbers of load op- 
posing components also increase the internal friction of the 
coupler. This lowers the ability of the coupler to efficiently 
track rescued kMT tips, as noted by the slower coupler veloc- 
ities in Fig. [6^. 

Overlap-opposing motors (such as, CENP-E kinesin) un- 
der tension load bring the passive binders at the tip of the 
kMT. We interpret these kinesins as components that increase 
the effective tension against the coupling which helps sustain 
tension-dependent suppression of depolymerization. Coupler 
velocities and Fb rea k do not show a significant difference as 
n + is varied under tension, Fig [6t. This suggests that these 
components do not significantly alter the nature of the kineto- 
chore attachment if sufficient numbers of passive binders are 
engaged. Kinesins in the our model serve to enhance the kMT 
tip tracking efficiency of the coupler, rather than significantly 
destabilize passive coupling. Finally, when plus end motors 
experience compressing forces, motor friction effects domi- 
nate coupler movement (SI). 

We conclude that the mean lifetime of the MT-kinetochore 
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attachment depends on a very delicate balance of forces and 
kinetic effects. Over a range of parameter values, our model 
reproduces the non-monotonic variation of the lifetime with 
load force that was observed in These attachment times 
give rise to optimal kMT- supported force ranges, indicating 
the existence of a force-dependent pathway for error cor- 
rection at kinetochores. We emphasize that the catch-bond- 
like phenomenon, that arises naturally from our microscopic 
model, is only one of the several distinct responses of the 
kMT-kinetochore coupler to load tension. We have also ex- 
plored other parameter regimes that correspond to (a) differ- 
ent lengths of the coupler formed by the binders, (b) differ- 
ent potential landscapes, and (c) different rates of growth and 
shrinking of microtubules. We have also investigated the ef- 
fects of a hybrid structure of the coupler that consists of an 
inner "active" layer of motor proteins and an outer "passive" 
layer of MT-binders. Depending on the parameter regime, the 
lifetime may appear to deviate from the qualitative features 
observed by Akiyoshi et al. lEHl . Some of the novel trends of 
variation observed in our analysis can be tested, in principle, 
by altering the size, composition, etc. of the single kineto- 
chore particle in in-vitro experiments. 
* Supplementary Information 
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SUPPLEMENTARY INFORMATION 

Kinetochore binder spacing: commensurate or incommensurate 
with MT lattice 

A key feature of our model is that the kinetochore fibrous 
corona is assumed to be densely populated by MT binder ele- 
ments. Even though these elements can be flexible, for a dense 
enough kinetochore region, we can assume a regularly spaced 
array of binder elements that are connected to a rigid coupler 
backbone. This approach allows us to construct an energy 
landscape as a function of the number of binders, Nb that can 
associate with the kMT. The binders are assumed to be equi- 
spaced and the spacing s between the successive binders is 
treated as a parameter that can attain arbitrary numerical val- 
ues. In this letter we have presented results only for the special 
case where s is commensurate with the MT lattice spacing, L 
However, in general, s need not be commensurate with MT 
lattice spacing 1291. One can interpret mconmesurate spacing 
as a more random arrangement of coupler heads. 

If the internal coupler friction increases significantly (mea- 
sured by b), the coupler loses its ability to adapt its position 
sufficiently fast with the changes in the position of the MT 
tip. That is why, in the main text, we refer to the strong fric- 
tion regime couplers as being practically static relative to the 
kMT lattice. Consequently, in this limit, detachment is possi- 
ble only because of the exit of the kMT tip from the coupler 
caused by its depolymerization. 

In the in-vitro experiments 112111 the binder heads might 
show heterogeneous binding energies, or even stochastic fluc- 
tuations in the numerical value of N^, which are ignored in 
our simple model. 

Microtubule Kinetic rates 



As in the previous work [I24L 12911 . we assume that the MT 
polymerization rate ot{x) inside the coupler (if non-zero) is 
small for all valuees of x, except for x = L where a = 
because of the lack of space between the MT tip and the kine- 
tochore wall. On the other hand, we assume that the depoly- 
merization rate f3 of the inserted MT slows down as tension 
force Fi oa d applied to the coupler increases. The polymeriza- 
tion and depolymerization rate functions used for the kMT tip 
kinetic rate r = a — /3(-Fi oa d) are given by 



a(x) 



1 + exp(-Ai(x - ai)) : 

/^load) = /W- AFload 



(7) 
(8) 



For a given F\ oa ^ the kMT can be in a state of either poly- 
merization (a > f3) or depolymerization (/? > a). For all the 
calculations shown here, a and f3 were such that the MT was 
in the state of depolymerisation when Fi oa d = 0. On gradu- 
ally increasing Fi oa d, the MT can switch to the state of poly- 
merisation because of tension-induced suppression of fi(F). 



In this model, we refrain from introducing stochastic transi- 
tions between states of polymerization and depolymerization 
using rescue and catastrophe frequencies. However, such an 
extension can be easily introduced in the model and we will 
explore it elsewhere. Of course, the introduction of stochas- 
tic catastrophe/rescue transitions in the kinetic kMT rates will 
add noise to the mean attachment time calculations. How- 
ever, as we'll report elsewhere, our key qualitative observa- 
tions made in this letter are not significantly affected by this 
noise. 



Load-Independent kMt rates for hybrid couplers. 



For this scenario we use 



fimax fin 



l + exp(Ai(x-/3i)) 



(9) 



in agreement with previous work 112911 . We note that in this 
case, r = a — f3 is adjusted such that the kMT experiences 
growth up to x < fii = 35 nm, with a > f3 m { n and subsequent 
destabilization for x > f3\ with a < /3 max . The rate a decays 
to zero if x > ai = L — I, where r transitions again, with 
f = —fimax- In the simulation results shown in Fig. S3 we use 
A™ = 27 s -1 , /3 max = 100 s~\ and a = 80 s -1 . 



Mean Attachment time calculations 



Here we give details of our calculations of the mean lifetime 
of the kMT-inetochore attachment within the framework of 
our model. We take a continuum approach JU]. Given that 
we start at a position x, the time required to exit from the 
boundary x = is denoted by T{x)\ it satisfies the delay 
differential equation 

-1 = j(-%(x) - F)d x T(x) + Dd 2 x T(x) + a(x)(T(x + £) - T{x)) 



(10) 



As shown in 12911 . in the slippery regime, the mean first exit 
time T{x) can be calculated using a reduced ordinary differ- 
ential equation, 



V\x, F, a, p)T x [x) + DT xx (x) = -1. 



where 



a(x) 



l2k B T 



(11) 



(12) 



V'(x,F,a,/3) = ^(-%(x)-F)-U^(f(x)/k B T)^(F)-a(x^ 

(13) 

and Iq (f(x) /ksT) is the integral form of the modified Bessel 
function of the first kind, which scales the effect of the MT 
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polymerization/depolymerization rates against the activation 
barrier heights in the potential energy well ^(x). The func- 
tion f(x), in the Bessel function, arises from a smoothed ap- 
proximation of the well where we use Sfr(x) = f(x)(l — 
cos(^p)) + h(x), with the linear terms f(x),h(x) chosen 
from a Fourier series fit to the potential energy well, (f(x) = 
jgx+l.b, h(x) = — g f). This approximation is made in order 
to simplify the numerical calculations. 

The boundary conditions are T(0) = 0, T'(L) = 0, where 
we have an absorbing boundary at x = and a reflecting 
boundary at x = L\ the reflecting boundary represents an im- 
penetrable physical barrier erected by the rigid kinetochore 
plate. For the mean attachment time, the Eq. SJTTT) yields the 
formal analytical solution 



In this case, 



T{x) 



1 

D 



dy 



exp(V(y)/D) 



[ exp(V(z)/D)dz. (14) 
Jy 



To get a closed form expression for the mean first pas- 
sage time we consider the special case of slippery limit where 
the activation barriers against transitions between successive 
binding sites are sufficiently low (b <C ksT) such that the po- 
tential well is well-approximated by Sk(x) = —ax, a = a/L 
Further, one can extend the binders to the coupler boundary, 
such that the well does not become flat inside the coupler (i.e., 
it remains linear). In this special case the contribution from 
the Bessel function is unity (2o(0) = 1), and the equation for 
the mean first passage time is simplified to 



'd-F 



^ ma xexp(-A(F)))r x (x) +DT xx {x) = -1, 

(15) 

where we have set a = without loss of generality. Solving 
this differential equation we get 



T(L) 



I? exp (—w) — 1 + w 



where 



w = 



L( -q + F + ^ max ^exp(-AF)) 



(16) 



(17) 



is a non-dimensional quantity. This formula provides a good 
estimate of the mean attachment times obtained from simula- 
tions in the limit b « ksT. 

In the strong friction regime, diffusion of the coupler along 
the kMT is insignificant and binding is very strong (a, b » 
kBT), so that movement of the kMT inside the coupler is only 
accomplished with the Poisson counting processes. One can 
alternatively think of this regime as a velcro or sticky regime, 
where the coupler is so strongly engaged with kMT binding 
sites that relocation relative to the kMT is not possible even 
under load, due to the high binding energy (coupler binding 
will eventually break under significant load, however for this 
limit coupler breaking would require load beyond the ranges 
tested in vitro). 



dx(t) ~ £dN r (t) 



(18) 



In the special case of a MT with a = 0, the overlap inside the 
coupler depends only on the depolymerization steps counted 
by the process dNp( F y For this process, the time interval 
between the occurrence of the consecutive depolymerization 
events are independent identically distributed random vari- 
ables with mean 1/(3(F). Therefore, the mean exit time for 
a fully engaged coupler with length L = N to d in the strongly 
bound regime is 



T(L) 



-AFioad 



(19) 



dadl-1 deletion in silico and force- velocity calculations 

Force dependent depolymerization is key in order to ob- 
serve the non-monotonic variation of the mean lifetime of the 
coupler with the load force. We test our model by reproduc- 
ing the results of the in-vitro experiments reported in 112 ill , 
where dadl — 1 components were deleted from the coupler. 
To mimic these deletions, we remove 5 binders and we set 
A = in order to remove lateral coupler cohesion provided by 
DAM ring components. In this scenario, the mean attachment 
times decrease monotonically with increasing load force (see 
Fig SO in full agreement with the corresponding observation 
of Akiyoshietal. ll2lll . 



x10 




N b =65,?i = 3 

N =60, X = 

- b 



200 




0.5 



1 1.5 

Force (pN) 



2.5 



FIG. 5. Mean attachment times in the slippery regime. We show 
calculations with Nb — 65 (solid line) and Nb — 60 (dashed line), 
using eq. fS[T6l). We also remove force dependent depolymerization 
in one case (A = 0), which completely removes range of forces for 
which attachment times increase. These results are used to reproduce 
Dadl-1 deletion studies from [21]. Figure inset. Plot of the kMT 
kinetic rates as a function of load force. Common parameters for all 
calculations shown are A = 3 pN -1 , /3 ma x = 140, a = 20, a = 0.4 
k£>T, k = 0.001a. 
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Velocity Reversals 

A suitably chosen value of the parameter N b or a can ensure 
that the mean attachment times exceed 1 min for the range of 
Fioad for which velocities are measured by simulations of eq. 
(1). In the inset of Fig. 3, we show coupler speeds for such 
a stably attached coupler in the slippery regime. The typical 
values of these speeds are comparable to the data of ref. [12111 . 

In the plot of the speed versus force, shown in Fig. 3inset, 
we see that the coupler experiences first a gradual slow-down 
and, beyond the velocity transition point, eventual speeding 
up as the tension load is increased further. The transition point 
corresponds toi ? c = l/ A = l/3 pN, which is consistent with 
the form of the /3(Fi oa d) function. We highlight that the value 
of F c can be increased only in the cases where there is suffi- 
cient a, N b to support attachment; in order for the nonlineari- 
ties to persist it is important that the depolymerization rate is 
slowed down within the range of forces that can be supported 
by the coupler. 



Numerical Simulations 

The parameters in the "intermediate regime" satisfy neither 
the conditions for slippery motion nor those for strong fric- 
tion. In this parameter regime the mean lifetime of the kMT- 
kinetochore attachments were computed by computer simula- 
tions of eq. (1) of the main text. The averaging was carried 
out over 500 trials for each fixed value of the force. Time step 
size had to be sufficiently small to ensure numerical stability 
and convergence. Exit time searches were extended up to 60 
min, which is well within the range of times explored experi- 
mentally in lEHl . In the stably attached cases, with high a, N b 
(or large attachment energy), numerical trials for exit times 
failed to yield an estimate, because the coupler remained en- 
gaged for times exceeding the allotted 60 min. Such an ar- 
tificial truncation of the lifetime distribution at 60 min gives 
rise to artificially lower estimates for the mean and standard 
deviation of attachment times at the exit time peak points in 
the intermediate regime; this is due to the smaller number of 
points used to gather exit time statistics. 



Active Coupler Interface 

Here we describe our approach for incorporating force- 
generation by active components at the coupler interface into 
the dynamical equations. We will refer to these elements as 
generic molecular motors, however, we note that in the case 
of kMT minus end motors (n_), then the force of the motors 
is equivalent to the force exerted by a ring that is being pushed 
by the curling plus end tips of kMTs. We note here that minus 
end motors will increase the kMT-coupler overlap, because 
these components walk toward the MT minus end, while the 
plus end of the kMT is the side that gains attachment with the 
coupler. On the other hand, for plus end motors, their action 



will decrease kMT-coupler overlap because the motors push 
toward the plus-end tip of the inserted kMT. 

For our purposes, we find it sufficient to adopt a stan- 
dard prescription for capturing the force-generation by the 
molecular motors that does not explicitly describe the detailed 
stochastic mechano-chemistry of the individual motors. 



Derivation of the model equations for the coupler composed of 
binders and motors. 



We start with force balance equation which does not include 
random fluctuations for the coupler overlap velocity 



dx(t) 
dt 



Vj 



MT 



%(x) - Fi 



load 



(20) 



where £ is the coupler effective drag coefficient and Vmt is 
the velocity of the kMT tip with respect to a space-fixed frame 
of reference. The active force term 



F A (x) = d m (x)(n-f- - 



(21) 



with the motor density function 



d m (x) = (x- N b s)(H(x - N b s) - H{x - N b s - L m )) 

(22) 

+ L m H(x — N b s — L m ), (23) 

where H(x) is the standard Heaviside step function and L m = 
8 nm corresponds to the total horizontal length of the cou- 
pler that can be populated by active components (in three- 
dimensions this corresponds to one layer of motors working 
around a kMT with 12 protofilament tracks, with one motor 
per track). 

As noted in the main text, linear force- velocity relations 
permit us to explicitly calculate the active component veloc- 
ity dependence for the total active components using /_ and 
/+. Note that minus here denotes minus end directed mo- 
tors/or protofilament curling that push to increase overlap and 
plus denotes plus end directed motors that work against cou- 
pler/microtubule overlap. For each case, we have 



f± 



1 



v± 



14 



(24) 



where F^ ax and V^ ax are the stall force and maximal veloc- 
ity for the plus-end directed and minus-end directed motors, 
respectively, whereas v± are the corresponding instantaneous 
velocities. Next we express v± in terms of dx/dt. Ignoring 
all the binder fibers, x = x t i P — x mo tor and hence 



dx 
~dt 



dx 



tip 



dt 
Vmt 



dt 



(25) 
(26) 



Note that, in the absence of plus-end-directed motors and ig- 
noring boundary conditions, the overlap can attain a stationary 
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value only if Vmt < 0; in this stationary state the depoly- 
merization of the tip would be balanced by the translocation 
driven by the minus-end-directed motors. Similarly, if only 
plus-end-directed motors are present, 



dx 
~dt 



Vmt - v + . 



(27) 



Substituting eqs. (l26t and (l27l) into eq. (f20l) we get 



dx 
~dt 



1 



-9'(x)-F ]oad + d(x){n-F r 



dx/dt — Vmt 



Vm 



x/dt^ 



Vj 



MT- 



(28) 



Regrouping the velocity terms we obtain the following 
equation for coupler overlap, 



1 



dx 

~db J(x) 
+ Vmt 



[-*'(x) - F load + d(x)(n-F max - n+F+ ax )] 

(29) 



where £(x) = £ + fi (x) 



- /i + (x) and 



± - d(x) 



n ±F^ ac 

^Vmax 



(30) 



end motors cause the minimum energy state to correspond to 
x = NfrS, which places the passive binders at the kMT tip. 
This effect causes the coupler to operate under a biased dif- 
fusion principle, especially when external loads are tension 
loads (Fioad > 0). 

The dynamics of the coupler described by eq. (l29t is fully 
deterministic. However, in reality, the kinetics of the kMT- 
kinetochore coupler are stochastic. Therefore, we now write 
down a stochastic differential equation (SDE) that would, 
upon averaging, correspond to the deterministic equations 
written above. Suppose over a small time interval St the num- 
ber of subunits (an a — f3 tubulin dimer) added and removed 
from the tip of the kMT by polymerization and depolymeriza- 
tion are dN r , an independent homogenous Poisson process. 

Moreover, we capture the effects of random Brownian 
forces through the noise W(t) which is assumed to be a Gaus- 
sian stochastic process. We distinguish this Gaussian process 
from the one used earlier in the absence of active force genera- 
tors. Since molecular motors are fueled by chemical reactions 
(e.g., ATP hydrolysis), this random force includes the effects 
of fluctuations both in the chemical reactions and mechanical 
stepping involved in each cycle of the individual motors. 

Thus, the equation for the coupler overlap reads 



dx 
~dt 



1 



[-*'(x) - Fioad + d(x)(n-F max - n + F+ ax j\ 

(31) 



CO*) 

■£dN r (t) + dW(t). 




x L 



FIG. 6. Potential binding energy modification upon active interface 
addition. 

We note that the potential well diagram in Fig. S0 demon- 
strates the effects of incorporating active components in the 
couplers. For minus-end motors the well is tilted in such a 
way so that the minimum energy state is located at x = L, 
for which the coupler is fully overlapped with the kMT. No- 
tice that the shape of the well indicates an increased internal 
friction arising not just due to the passive binder heads (lo- 
cal wells) but also due to the concavity of the modified well 
due to the minus end motors. On the other hand, for plus 



Force velocity relations for a hybrid coupler 

Here we show force- velocity calculations for a coupler mo- 
tor that is engaged with varying densities of plus-end directed 
or opposing motors and the inserted kMT kinetic rates do not 
have load dependence. 




-5 -4 -3 -2 -1 
Force (pN) 



FIG. 7. Average speed versus force plot for a stably bound coupler 
with varying densities of plus-end directed motors. Inset. Stall force 
calculations for each motor density inside the coupler. The param- 
eters are N b = 52, A = 3, a = k B T, k 0.01a, /3 min 27 s~\ 
(3 max 100 s" 1 , and a = 80 s" 1 . 



Force velocity relations for a hybrid coupler 

In order to challenge active components that work to de- 
crease the overlap (such as CENP-E kinesin), in Fig. S|7]we 
tested the model with varying n + and Fi oad < 0. In these 
cases, Fi oad < is not a tension load, so we omit kMT f3 
dependence on Fi oa d. As shown in Fig. $71 when the ac- 
tive motor component pulls to oppose coupler overlap against 
compressing load, the internal friction is directly increased by 
the action of the motor, resulting in an overall loss of cou- 
pler tracking ability (however attachment is maintained, the 
motors become stalled), as shown by decreasing velocities in 

Fig. sm 

ERROR CORRECTION AT KINETOCHORES AND 
CATCH-BOND-TYPE MECHANISMS 

The non-monotonic attachment times which arise for vari- 
ous parameter ranges for this coupler model indicate that there 
are force ranges for which longer attachment times are fa- 
vored. We interpret the regions with peak attachment times 
as stable-attachment force regions. Thus the force ranges that 
give long T(L) can be thought of as optimal force ranges for 
which kinetochores support attachment. This force- selective 
mechanism is particularly important in the context of error 
correction of kMT attachments at kinetochores; many erro- 
neous kMT attachments may not provide sufficient tension 
force on a kinetochore due to the geometry of the connec- 
tion. A purely force-mediated error-correction mechanism at 
kinetochores has important implications because it expands 
the role of kinetochore couplers beyond generation of move- 
ment, to also checking the quality of the kMT attachment. 

There are various mechanisms for catch-bonds that explain 
how such a bond can become stronger under force (for review 
see [? ]). In the general context of MT-kinetochore attach- 
ments, references to 'catch-bonds' are made to account for 
the observed increase in the attachment life time when sub- 
jected to tensile force. We have offered a conceptual scenario 
where the stronger is the tension the smaller is the curvature 
of the splaying kMT tips and, as a consequence, the slower is 
the depolymerization of the kMT and, hence, the longer is the 
life time. 

To implement this conceptual model quantitatively, we have 
assumed that the depolymerization rate is proportional to the 
free depolymerization rate (or attempt frequency) and expo- 
nentially related to the height of the energy barrier along the 
unbinding pathway. We have decoupled the coupler binding 
energy by keeping the energy function *ff(x) and modified the 
depolymerization rate by a rescaled load force term, Fi oa dA. 
The two can be combined if we use a dynamic potential en- 
ergy landscape which takes into account protofilament curling 
energies (\P (#,£)), however, we claim that the end effect is 
similar to the simple model we present here. 
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PARAMETER VALUES 



Parameter Description 


Symbol 


Values Tested 


MT binding site spacing 


£ 


8/13 nm [23J 


Maximal coupler length 


L 


50 nm [? ] 


Maximal number of coupler binders 


N b 


15-65 


Maximal Depolymerization Rate 


fimax 


100 - 350 s' 1 [23J 


Maximal Polymerization Rate 


a 


20 - 50 s -1 [24] 


MT lattice/binder binding energy 


a 


0Ak B T - 3k B T 


Binder activation barrier 


b 


0.001a -0.4a 


Critical Depolymerization Force 


F c = 1/A 


0.3 - 5pN 


Polymerization Decay Position 


an 


L - £ nm [29J 


kMT-rate transition steepness 


Ai 


100 nm" 1 [29] 



